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FOREWORD 


This manual, which doc\aments the THERMAL RESPONSE FAC- 
TOR ROUTINE (RESFAC - version 2) used, by NASA'S ENERGY COST 
ANALYSIS PROGRAM (NECAP) , is a suppliment to the NECAP ENGI- 
NEERING MANUAL (TM 83240) and the NECAP ENGINEERING FLOWCHARTS 
MANUAL (TM 83242) . 

The calculation sequences and flowcharts are given in 
the same format that is used by the NECAP ENGINEERING MANUAL 
and the NECAP ENGINEERING FLOWCHARTS MANUAL . 

This version of RESFAC is the result of modifications 
and enhancements made to the original RESFAC, Although the 
RESFAC programming code has changed, the way in which the 
routine is accessed and used by NECAP is still the same. 
Therefore, the RESFAC usage information presented in the 
NECAP INPUT MANUAL (TM 83239) and the NECAP USERS MANUAL 
(TM 83238) is still accurate. 
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OBJECTIVE AND DESCRIPTION 


The Response Factor Program (RESFAC2) generates the 
time series of heat transfer factors (called response fac' 
tors) required to accurately determine the transient flow 
of heat into, through, and out of interior/exterior buil- 
ding surfaces as they react to temperature differences 
across them. 

These response factors are a function of the type 
of materials used and their order of placement in the 
surface. It is required that the following material 
properties be known for each material layer: 
for non-air layer- 


1. 

XL, 

thickness (FT) 

2. 

XK, 

thermal conductivity (BTU/HR-FT-“F) 

3. 

D , 

density (LB/FT3) 

4. 

SH, 

specific heat (BTU/LB-®F) 


for air layer- 

5. RES, resistivity (HR-Ft2-°f/BTU) 

The sequencing of each layer's material character- 
istics is important. It must follow the way the surface 
is layered from the outside material to the inside. 
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DEFINING THERT4AL RESPONSE FACTORS 


Thermal response factors can be defined as time depen- 
dent heat transfer coefficients which depict an object's 
action upon a particular outside surface/inside surface 
temperature difference over a period of time. This action 
is composed of three terms — the X, Y, and Z response factors. 

Using the heat £l\ix equation Q « U along with the 

response factors, the time dependent heat flux values for 
the planar construction shown in Figure 1 can be obtained. 



INSIDE SURFACE 

QI 2 
(BTU/HR-FT ) 

TI 

(“F) 


FIGURE 1 


The outside surface heat flux (QO) at any time t is 
calculated by : 


NRT NRT 

°0t - 

i=0 i=0 

where NRT is the number of response factor terms calculated, 

TO, . and TI. . are outside surface and inside surface tem- 
t-j. t-i 

peratures at times t-i hour. 
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The inside surface heat flux (QI) at any time t is cal- 
culated by ; 

NRT NRT 

Qi^ = 

i=0 i=0 

In order for the steady state heat conduction situation 
QO^ = QI^ = U(TO - TI) 

TO = TO^ for any time t 

TI = TI^ for any time t 

to be satisfied, the individual summation of the X, Y, and Z 
response factors must equal the overall heat transfer coef- 
ficient U : 


u = 

OB 

Ex, 

s: 

OB 

Ex, 


OB 

Ex, 


i=l 


i=l 


i=l 


3 



ALGORITHMS 


The new Response Factor Program is comprised of five 
program modules : - 

Main program - RESFAC 

Subroutines - RESIDU, ROOTS, and THREBS 
Function - BMEQ 

The following pages present input/output, calculation 
sequences, and flowcharts for each of these modules. 
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RESFAC 


A program which determines the thermal response factors 
for a particular roof, wall, or floor construction. 


INPUT 


N 

XL 

XK 


SH 


HES 


DEFC 


nvunber of material layers 

thickness of each layer (FT) 
if air layer then XL >= o.O 

c 

thermal conductivity of each layer (BTU/HR-FT-F) 
if air layer then b o.O 

density of each layer (LB/FT*) 
if air layer then D = o.O 

o 

specific heat of each layer BTU/Lb-F) 
if air layer then SH = 0.0 

thermal resistance of each AIR layer 
(HR-FT*-*t/BT0) 

if non-air layer then RES = 0.0 

alpha-numeric description of each layer, 
maximum of 30 characters 


OUTPUT 


RFX ^ 
RFY > 
RFZ J 

NRFT 

Rl 


Response factors series 

number of response factor terms 
common ratio 


CALCULATION SEQUENCE 

1. Read number of material layers M 

2. For i = 1, M 

a) read and echo 
1 ) 

2 ) conductivity XK^ 

3 ) density 

4 ) specific heat SH^ 


thickness XL^ 
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5 } resistemce RES^ 



6) 

alpha-numeric description DEFC^ 

b) 

for 

each non-air layer, calculate: 


1) 

its resistance R^ *= 


2.) 

BETA^ = XL^*SQRT(D^*SHj^/XK^) 

c) 

for 

each air layer: 


1) 

set Rj^ •= RES^ 


2) 

BETA^ « 0.0 


3. Call subrou-txne EESIDU 1:o compu1:e the thermal conductance 
RKO emd the derivatives o£ matrix elements A, B, and D 
(variables RKl, RHl, and HIi4) for the condition where P, 
of B(P)r approaches 0.0 

4. Call subroutine ROOTS to calculate the roots of B(P) = 
0.0 and the number of roots NRT 

5. For i = 1, NRT 


1 ) call subroutine THREBS to calculate matrix 
elements A, D and. the derivative of matrix 
element B (variables Bl, B2, and BP3) 

2) set RKK(i,2) = l-0/BP3/ROOT(i)^ 

RKK(i,l) = RKK(i,2)*Bl 
RKK(i,3) = RKK^,2)*B2 

6. For i=l,100: 

a) for j = 1, NRT: 

if (ROOT(j) > 30.0) go to 6b 
BETAY = EXP(-ROOT(j) *i) 

A = A + RKK(j,l)*BETAY 
B = B + RKK(j ,2) *BETAY 
C = C + RKK(j,3) *BETAY 

b) A = A + RKl + (RM4*RK0) + (i*RK0) 

B = B + (i*.RK0) +. RKl 

C = C + RKl + (RM1*RK0) + (i*RK0) 

c) if Xi = 1) : 

AA = A 

BB = B 

cc = c 
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d) if .(i = 2) : 

AA *= A-2.0*X 

BB = B-2.0*y 

CC = C-2.0*Z 

e) if ( > 2) : 

AA = A-2.0*X + FA 

BB = B-2.0*y + FB 

CC = C-2.0*Z + FC 

'£) set and echo response factors: 

RFX(i) = AA 
RFy(i) = BB 
RFZ(i) = CC 

g) if (i < 3) go to 6k 

h) if CABS(XX/FAA-AA/30C) > 1.0E“5) go to 6k 


i) if jCABS(yy/FBB-BD/yy) > 1.0E-5) go to 6k 


j) if (ABSCZZ/BCC-CC/ZZ) < l.OE-5) go to 7 


k) if (ABS(AA) < l.OE-8) go to 7 

J^) set: 

FA = X 

FB = y 

FC = Z 
FAA = XX 

FBB = yy 

FCC = ZZ 
XX = AA 
yy = BB 
ZZ = CC 
X = A 

y = B 

z = c 

go to 6 

7. Set and echo common ratio: 

Rl = EXP (-ROOT (1)) 

8. End 
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RESIDU 


A subroutine which determines thermal conductance and the 
derivatives of matrix elements A, B, and D for the condition 
where the root W approaches 0.0, 


INPUT 


RR 

BETA 

M 

OUTPUT 

RKO 

RKl 

RMl 

RM4 


each layers resistance 
each layers BETA 
number of layers 

thermal conductance 
derivative of matrix element B 
derivative of matrix element A 
derivative of matrix element D 


CALCULATION SEQUENCE 


1. Sum all layer resistances and compute thermal 
conductance RKO = 1.0/R total 

2. For i = 1,M 

a) P = BETA(l) *BETA(1) 

R = RR(1) 

b) if (i > 1) : 

A = 1.0 

B = R 
C = 0.0 
D = 1.0 
go to 2e 

c) if (i = 1) i 
A = P/2,0 

B = R*P/6.0 
C = P/R 
D = P/2.0 

d) if (M =1).; 

RMl = A 

RKl = -RK0*RK0*B 
RM4 = D 
go to 4 
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e) for j = 2 , M 

1) P = BETA( j)*BETA( j) 
R = RR (j ) _ 

if i = j go to 2e2 

AA = 1.0 

BB = R 

CC = 0.0 

DD = 1.0 

go to 2f 

2) AA = P/2.0 
BB = P*R/6.0 
CC = P/R 

DD = P/2.0 

f) FI = A*AA+B*CC 
F2 = A*BB+B*DD 
F3 = C*AA+D*CC 
F4 = C*BB+D*DD 
A = FI 

B = F2 
C = F3 
D = F4 
go to 2e 

g) RMl = RMl+Fl 
RM2 = RM2+F2 
RM4 = RM4+F4 
go to 2 

3. RKl = -RK0*RK0*RM2 

4. End 
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ROOTS 


A siibroutine which determines tlie roots of matrix element 
B(P) = 0.0 


INPUT 

BTA 

R5 

N 

OUTPUT 

ROT 

L 


each layers BETA 

each layers resistance 

number of material layers 

root array 
number of roots 


CALCULATION SEQTJENCE 

1. Set PI = 3.141592654 

2, For i = 1, N: 


a) if BTA(i) >0.0 then calculate the zero points 
of this layers sine euid cosine components: 

1) set R = 1.0 

2) for j = 1, 99, 2 

BM(j, i) = (R - 0.5)*PI/BTA(i) 

BM(j + 1, i) = R*PI/BTA(i) 

R = R + 1.0 

3. For values >0.0, sort the BM arrays' values into the 
possible root interval array PINT in ascending order. 

4. For i = 1,1000 determine the roots such that no 
root >30.0: 

a) set G1 = PINT(i) and calculate Fl(Gl) using 
f\inction BMEQ 

1) if (ABS(Fl) < l.OE-8) go to 4e 

b) set G2 = PINT(i + 1) and calculate 
F2 (G2) using function BMEQ 

1) if (ABS(F2) < l.OE-8) go to 4e 

c) set the interval length D = G2 - G1 
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d) sectionaiize each interval into 1/lOth segments 
and check each segment to determine i£ there is 
a root there; 

1) FK = FK + 1.0 

if (FK > 10.0) go to 4 

2) set DD = PINT(i) + (D*FK/10.0) and 
calculate FD(DD) using function BMEQ 

3) SM = FP*F1 

if (SM <= 0.0) go to 4d5 

4) G1 = DD 
FI = FD 
go to 4dl 

5) GX = G1 

FX = FI 

G2 = DD 

F2 = FD 

go to 4f 

e) L = L + 1 

1) if (ABS(Fl) <=1.0E-8)ROT(L) = Gl**2 

2) if (ABS(F2) <= l.OE-8 and RqT(L) = 0.0) 

ROT(L) = G2**2 

3) if (RCr(L) > 30.0) go to 5 
else go to 4 

f) use secant method to find root 

1) G3 = G2 - (F2*(G2-G1)/(F2-F1) ) 

2) if (G3 < GX or G3 > DD) then secant method fails 
so use bisection method to locate root - go to g 

3) if fcl3 >GX and G3 <DeJ then calculate FRT(G3) 
using function BMEQ 

4) if (ABS(FRT) <= l.OE-8) go to 4j 

5) G1 = G2 

G2 = G3 

FI = F2 

F2 = FRT 

go to 4fl 
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g) 

set G4 = 0.5* (G1 + G2) 
F4(G4) using function 

and calculate 
BMEQ 

h) 

if (ABS(F4) <= l.OE-8) 

G3 = G4 


1) go to 4j 


i) 

F5 = F1*F4 



if (F5 < 0.0) G2 = G4 

if (F5 > 0.0) G1 = G4 and Fl = F4 

go to 4g 

j) L = L + 1 

. ROT (I) = G3**2 

if (ROT(L) > 30.0) go to 5 
G1 = DD 
Fl = FD 
go to 4dl 

5) End 
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THREBS 

A siibroutine which determines matrix elements A, D,and the 
derivative of matrix element B at the root W. 


INPUT 

RR 

BETA 

W 

M 

OUTPUT 

B1 

B2 

BP3 


each layers resistance 
each. layers BETA 
root 

momber of layers 

matrix element A 
matrix element D 
derivative of matrix element B 


CALCUALTION SEQUENCE 
1. For i = 1, M 

a) P = SQRT(W) *BETA(i) 

R = RR(i) 

ALPHA = BETA(i) 

SQ = SQRT(W) 

b) if (P = 0.0) ; 

AL(i) = 1.0 
BL(i) = R 
CL(i) = 0.0 
DAL(i) = 0.0 

DBL(i) = 0.0 

DCL(i) = 0.0 

C) if (P > 0. 0) : 

AL(i) = COS(P) 

BL(i) = R/P*SIN(P) 

CL(i) = -P/R*SIN(P) 

DAL(i) = ALPHA*SIN(ALPHA*SQ)/2. 0/SQ 
DBL(i) = -R*COS (ALPHA*SQ)/2. 0/W+R*SIN (ALPHA*SQ)/ 
ALPHA/2 . 0/SQ/SQ/SQ 

DCL(i) = ALPHA*ALPHA*COS(ALPHA*SQ)/2.0/R+SIN 
(ALPHA*SQ) /2 . 0/SQ*ALPHA/R 

d) go to 1 
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2. If (M = 1) : 

BP3 = DBL(l) 

B1 = DAL(l) 

B2 = AL(1) 
go to 7 

3. For i = 1, M 

a) DA = 1.0 
DB = 0.0 
DTA =0.0 
DTB =0.0 

b) for j = 1, M 

1) if (i = j) s 

DTA = (DA*DAL(j)) + (DB*DCL(j)) 
DTB = (DA*DBL(j)) + (DB*DAL(j)) 
go to 3b 3 

2) DTA = (DA*AL(j)) + (DB*CL(j)) 
DTB = (DA*BL(j)) + (DB*AL(j)) 

3) DA = DTA 
DB = DTB 
go to 3b 

c) BP3 = BP3 + DTB 

d) go to 3 

4. A = AL(1) 

B = BL(1) 

C = CL(1) 

D = AL(1) 

5. For i = 2, M 

a) TA = (A*AL(i)) + (B*CL(i)) 

TB = (A*BL(i)) + (B*AL(i)) 

TC = (C*AL(i)) + (D*CL(i)) 

TD = (C*BL(i)) + (D*AL(i)) 

A = TA 
B = TB 
C = TC 
D = TD 

b) go to 5 

6 . B1 = TD 
B2 = TA 

7. End 
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BMEQ 


A function which determines the matrix element B at the 
root guess W. 

INPUT 

W : root guess 

BETA : each layers BETA 

R : each layers resistance 

M : number of layers 

OUTPUT 

BMEQ : matrix element B 

CALCULATION SEQUENCE 

1. For i = 1, M 

A) if (BETA(i) = 0.0) : 

AL(i) = 1.0 
BL(i) = R(i) 

CL(i) = 0.0 
to to 1 

b) P = W*BETA(i) 

AL(i) = COS(P) 

SN = SIN(P) 

BL(i) = R(i)*SN/P 
go to 1 

2. If (M = 1) : 

TB = BL(1) 
go to 5 

3. A = AL(1) 

B = BL(1) 

4. For i = 2, M 

TA = (A*AL(i)) + (B*CL(i)) 

TB = (A*BL(i)) + (B*AL(i)) 

A = TA 
B = TB 


5. BMEQ = TB 

6 . End 
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